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Abstract. We give sufficient conditions under which the convergence 
of finite difference approximations in the space variable of possibly de- 
generate second order parabolic and elliptic equations can be accelerated 
to any given order of convergence by Richardson's method. 



1. Introduction 



This is the third article of a series studying a class of finite difference 
equations, related to finite difference approximations in the space variable 
of second order parabolic and elliptic PDEs in M. d . These PDEs are given 
on the whole M. d in the space variable, and may degenerate and become first 
order PDEs. Denote by Uh the solutions of the finite difference equations 
corresponding to a given grid with mesh-size h. By shifting the grid so that 
x becomes a grid point we define Uh for all x £ M d rather than only at the 
points of the original grid. In [5] and [6], the first and second articles of the 
series, we focus on the smoothness in z of rather than their convergence 
for h — ► 0. The main results in [5] and [6] give estimates, independent of 
h, for the first order derivatives Duh and for derivatives D k Uh in x of any 
order k, respectively. 

In the present paper one of our main concerns is the smoothness of the 
approximations Uh in (x,h). In particular, we are interested in the conver- 
gence of Uh, and their derivatives in x, in the supremum norm, as h — > 0. We 
give conditions ensuring that for any given integer k > the approximations 
Uh admit power series expansions up to order k + 1 in h near like 



and such that the coefficients are bounded functions of (t, x) € [0, T] x M. d 
for fixed T > in the case of parabolic equations, and, with the exception 
of rh, are independent of h. This is Theorem 12.31 our first result on Taylor's 
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formula for Uh in h. We obtain it by proving first Theorems 12.11 and 12.21 
below on the solvability of the PDE that is being approximated, and of a 
system of degenerate parabolic PDEs, respectively, for the coefficients 
j = 0, . . . , k. Of course, is the true solution of the corresponding PDE. 
The remainder term rt satisfies a finite difference equation, with the same 
difference operator appearing in the equation for Uh, and we estimate 
by making use of the maximum principle enjoyed by this operator. This 
is a standard approach to get power series expansions for finite difference 
approximations in general, and it works well in many situations, when suit- 
able results regarding the equations for the coefficients are available. 
In our situation it requires some facts either from the theory of diffusion 
processes or from the theory of degenerate parabolic equations. However, 
we do not use any facts from these theories. We prove Theorem 12.11 and 
hence Theorem 12. 2\ relying on results on finite difference schemes, obtained 
in [6] by elementary techniques. It is worth saying that since long ago finite 
difference equations were already used to prove the solvability of partial dif- 
ferential equations (see, for instance, [8] and [9]). Our contribution lies in 
considering degenerate equations. 

After establishing the expansions of in h not only we obtain the pos- 
sibility to prove the convergence of Uh to the true solution in the sup norm 
as h — > but also the possibility to accelerate it to any order under appro- 
priate assumptions. We prove the latter by taking linear combinations of 
finite difference approximations corresponding to different mesh-sizes. This 
method is especially effective when many of the coefficients in the expansion 
of Uh are zero. These results are given by Theorem 12.211 and Corollary 12.81 
Their counterparts in the elliptic case are presented by Corollary 13.71 

The idea of accelerating the convergence of finite difference approxima- 
tions in the above way is well-known in numerical analysis. It is due to 
L.F. Richardson, who showed that it works in some cases and demonstrated 
its usefulness in [15] and [16]. This method is often called Richardson's 
method or extrapolation to the limit, and is applied to various types of ap- 
proximations. The reader is referred to the survey papers [2] and [I] for a 
review on the history of the method and on the scope of its applicability 
and to the textbooks (for instance, [10] and [TT]) concerning finite difference 
methods and their accelerations. 

We are interested in approximating in the sup norm not only the true 
solution but also its derivatives. Note that even if the coefficients u^' are 
bounded smooth functions of (t, x), the derivatives D k Uh of Uh in x may not 
admit similar expansions, since the derivatives of may not be bounded 
in h near 0. Note also that the bounds on the sups of u^' and generally 
depend on T, and may grow exponentially in T. This becomes a big obstacle 
on the way of extending our results to the elliptic case. 

Our next result on power series expansions, Theorem 12.71 improves the 
previous theorem in two directions. It gives sufficient conditions such that 
for any given integer k > 



ACCELERATED SCHEMES 



3 



(a) D k Uh admits an expansion similar to (jl.ll) . 

(b) the bounds on the coefficients are independent of T. 

Having (a) we can approximate the A;-th derivatives of the true solution 
by D k Uh with rate of order h and accelerate the rate under appropriate 
assumptions. We can also approximate the £>th derivatives of the true 
solution with finite difference operators in place of D k applied to Uh, which 
is more convenient in applications because it does not require computing 
the derivatives of Uh- 

We ensure (a) and (b) by relying heavily on derivative estimates, indepen- 
dent of T, obtained in [5] and [6] for solutions of finite difference equations. 
Property (b) of the expansions allows us to extend Theorem 12.71 to the el- 
liptic case. This extension is Theorem 13.51 

As a consequence of the derivative estimates proved in [6] we obtain also, 
see Theorem 12.91 below, estimates, independent of h and T, for the deriva- 
tives of Uh in x and h. Clearly, Theorem 12.91 immediately implies Taylor's 
formula for in h, up to appropriate order, with bounded coefficients. It 
is interesting to notice that the converse implication does not hold: If for 
k > 1 the function Uh admits a power series expansion up to order k + 1 in 
h near with bounded coefficients, it does not imply, in general, that the 
derivative of Uh in h up to order k + 1 are bounded functions. That is why 
Theorem 12.71 does not imply Theorem 12.91 an d the latter implies the former 
only if condition (i) in Theorem 12.71 is satisfied. Additional information on 
the behaviour of the derivatives of Uh in x and h when h is near is given by 
Theorem 12. Hi The corresponding result in the elliptic case is Theorem 13.41 

In this article we are working with equations in the whole space having 
in mind considering equations in bounded smooth domains in a subsequent 
article. Still it may be worth noting that the results of this article are 
applicable to the one dimensional ODE 



The point is that one need not prescribe any boundary value of u at the 
points ±1 and if one considers this equation on all of K, the values of its 
coefficients and / outside (—1, 1) do not affect the values of u(x) for \x\ < 1. 

2. Formulation of the main results for parabolic equations 

We fix some numbers ho,T € (0, oo) and for each number h 6 (0, ho] we 
consider the integral equation 



for u, where g(x) and f(s, x) are given real-valued Borel functions of x € M. d 
and (s,x) G Ht = [0, T\ x M. d , respectively, and Lh is a linear operator 
defined by 



(1 — x 2 ) 2 u"(x) — c(x)u(x) = f(x) 



xe (-1,1). 




(2.1) 



L h ip(t, x) = L h tp(t, x) - c(t, x)ip{x) 



(2.2) 
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L° h p(t,x) = - ^ q\(t,x)5 h) \<p(x) + ^ pa^Mm^)' (2.3) 



AgAi AgAi 



for functions <p on R . Here Ai is a finite subset of M. such that ^ Ax, 

5h,\p{x) = -Affix + hX) - p(x)), AgAi, 

Qx(t,x) > 0, p\(t,x), and c(i, x) are given real-valued Borel functions of 
(t,x) eH QO = [0,oo) x R d for each A G A x . Set |Ai| 2 = £ AgAi |A| 2 . 
As usual, we denote 

D a = D?...D% 1 , Di = |o|=^ai, =DiDj 

% i 

for multi-indices a = (a±, . . . ad), ati G {0, 1, . . . }. For smooth tp and integers 
k > we introduce D k ip as the collection of partial derivatives of tp of order 
k, and define 

\D k p\ 2 = \D a p\ 2 , [p]k = sup \D k <p(x)\, \<p\ k = Y^[<P]i- 

\a\=k ^ i<k 

For functions iph depending on h G (0, Hq\ the notation D^tph means the k- 
th derivative of tp in h. For Borel measurable bounded functions tp = tp(t, x) 
on Ht we write ip G 25 m = SB™ if, for each t G [0, T], tp(t,x) is continuous 
in M. d and for all multi- indices a with \a\ < m the generalized functions 
D a ip{t,x) are bounded on Ht- In this case we use the notation 



sup V \D a rp(t,x) 



|a|<m 

This notation will be also used for functions ip independent of t. 
Let m > be a fixed integer. We make the following assumptions. 

Assumption 2.1. For any A G Ai, we have p\,q\,c, f,g G s B m and, for 
k = 0, m and some constants M k we have 



sup { £ (\D k qx\ 2 + \D k p x \ 2 ) + \D h c\ 2 ) < M 2 k . (2.4) 

Ht AGA! 



Remark 2.1. By Theorem 2.3 of [5] under Assumption 12.11 for each h G 
(0, ho], there exists a unique bounded solution of (|2.ip . this solution is 
continuous in Ht, and all its derivatives in x up to order m are bounded. 
Actually, in Theorem 2.3 of [5] it is required that the derivatives of the data 
up to order m be continuous in Ht, but its proof can be easily adjusted to 
include our case (see Remark 12.61 below) . 

Naturally, we view (12. ip as a finite difference schemes for the problem 
d 

—u(t,x) = £u(t,x) + f(t,x), t G (0,T], x G R d , (2.5) 
u(0,x)=g(x), x G M d , (2.6) 
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where 

d d 

£:=1E E ri\\;\jl),l)j + Y,P^ XiDi ~ C - ( 2 - ? ) 

AGAi i,j=l AgAi i=l 

By a solution of (|2,5p - (l2.6p we mean a bounded continuous function u(t, x) 
on Ht, such that it belongs to 25 2 and satisfies 

u{t,x) = g{x) + f [Cu(s,x) + f(s,x))ds (2.8) 
Jo 

in Ht in the sense of generalized functions, that is, for any t £ [0, T] and 

"»OC 



(j)(x)u(t,x) dx = / <f)(x)g(x) dx + / / (f>(—cu + f)(s,x)dxds 

JR d Jo JR d 

i-t /> d d 

+ / / § ^ (l ^ qxXiXjDjDjU + y^ j px\ i D i u)(s,x) dxds. (2.9) 

Observe that if u € OS 2 , then (|2.9|) implies that (|2.8p holds almost everywhere 
with respect to x and if u € 2$ 3 then the second derivatives of u in x are 
continuous in x and (|2.8p holds everywhere. 

The reader can find in [7] a discussion showing that in all practically 
interesting cases of parabolic equations like (|2.8j) the operator £ can be 
represented as in (|2.7p . so that considering operators l? h in form (|2.3|) is 
rather realistic. 

The following theorem on existence and uniqueness of solutions is a clas- 
sical result (see, for instance, |12j . |13| . [14] ) which we are going to obtain 
by using finite-difference approximations. 

Theorem 2.1. Let A ssumption \2.1\ hold with m > 2. Then equation (|2.8p 
has a unique solution € OS 2 = *B T . Moreover, € 03™ and 

\\u^\\ m <N{\\f\\ m + \\g\\ m ), (2.10) 

where N is a constant, depending only on d, m, |Ai|, Mo,. . . , M m , and T. 

Observe that this result is rather sharp in what concerns the smoothness 
of solutions, which is seen if all the coefficients of C are identically zero and 
/ is independent of t in which case the solution is tf(x) + g(x). 

The existence part in Theorem 12.11 is proved in Section [6] and uniqueness 
in Section HI 

In Section [6] a repeated application of this theorem allows us to prove a 
result on the solvability of (I2.13|) below. First introduce 

^ ■= (iTDW E ^ + *+T E (2-11) 

where 



(t+l)(i+2) 

AgAi AgAi 



9 W := E ^i D W ( 2 - 12 ) 
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is the derivative of <p in the direction of A. Consider the system of equations 

ft 3 

u®(t,x)= / (Cu^\s,x) + y^C]C^u^- i \s,x))ds, (2.13) 
Jo i=i 



(t,x) G H T , j = l,...,k. 



Remark 2.2. Quite often in the article we use the following symmetry con- 
dition: 

(S) Ai = — Ai and q\ = q-\ for all A 6 Ai. 
Notice that, if condition (S) holds, then 

h~ l ^2 q\(t,x)5 h ,x<p(x) = (1/2) ^2 q\(t,x)A h , x <p{x), 
AeAi AeAi 

where 

Ah t x<p(x) = h~ 2 ((p(x + hX) — 2ip(x) + (p(x — h\)). 

Theorem 2.2. Let k > 1 be an integer, (i) If Assumption POl is satisfied 
with m > 3k + 2, then (|2.13p has a unique solution {u^'}j =1 , such that 

U U) € gjm-s^ \\ u U)\\ m _ 3 . < TVdl/IU + \\g\\ m ) (2.14) 

for j = 1, . . . ,k. 

(ii) If the symmetry condition (S) holds and Assumption \2.1\ is satisfied 
with m > 2k + 2, then (12.13P has a unique solution {u^}j =1> such that 

U U) € fcm-iy h U)\\ m _ 2j < N{mm + || ff || m) (2 . 15) 

for j = 1 . . . , k. In addition, if 

P-\ = ~PX, for AG Ai, (2.16) 

then 

u ij) = 0, (2.17) 

for odd numbers j < k. 

In all cases the constants N depends only on d, m, |Ai|, M$, . . . , M m , 
and T. 

The next series of results is related to the possibility of expansion 

1 j 

u h (t,x) = u (0) (t,x) + V —u^(t,x) + h k+1 r h (t,x), (2.18) 

for all (t, x) G Ht and h G (0, ho], where Uh is the unique bounded solution of 
(|2.ip (sec Remark [2. ip and Th is a function on Ht defined for each h G (0, ho] 
such that 

\r h (t,x)\ < N{\\f\\ m + \\g\\ m ) (2.19) 

for all (t,x) G Ft, h G (0, /i ]. 
Introduce 

Xh,x = q\ + hp x . 
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Assumption 2.2. For all (t, x) G Hr, h G (0, /ioL and A G Ai, 

X M (t,x)>0. (2.20) 

Assumption 2.3. We have 

^2\q x (t,x) = for all (t,x) G H T . 
AGAi 

Notice that condition (S) is stronger than Assumption 12.31 

Theorem 2.3. Let Assumption \2.1\ with rri > 3 and Assumption \2.2\ hold. 
Let k > be an integer. Then expansion (|2.18p ZioWs wii/i satisfying 
(|2. 19[) . provided one of the following conditions is met: 

(i) m >3k + 3 and Assumption \2.3\ holds; 

(ii) m > 2k + 3 and condition (S) holds; 

(Hi) k is odd, m > 2k + 2, and conditions (S) and (|2.16p are satisfied. 
In each of the cases (i)-(iii) the constant N depends only on d, m, |Ai|, 
Mq, . . . , M m , and T. In case (Hi) we have = for all odd j in expansion 

We prove this theorem in Section [3 The following corollary is one of the 
results of [3] proved there by using the theory of diffusion processes. We 
obtain it immediately from case (iii) with k = 1. Of course, the result is 
well known for uniformly nondegenerate equations but we do not assume 
any nondegeneracy of C, which becomes just a zero operator at those points 
where q\ = p\ = c = 0. 

Corollary 2.4. Let conditions (S) and (|2.16p be satisfied. Let Assumption 
\2.1\ with m = 4 and A s sumption 1 2. 2\ hold. Then we have \uh — uq\ < Nh? . 

Actually, in [3] a full discretization in time and space is considered for 
parabolic equations, so that, formally, Corollary 12.41 does not yield the cor- 
responding result of [3]. On the other hand, a similar corollary can be 
derived from Theorem 13.51 below which treats elliptic equations and it does 
imply the corresponding result of [3]. It also generalizes it because in [3] one 
of the assumptions, unavoidable for the methods used there, is that q\ = r\ 
with functions r\ that have four bounded derivatives in x, which may easily 
be not the case under the assumptions of Theorem 13.51 

To formulate our main result about acceleration for parabolic equations 
we fix an integer k > and set 

k 

u h = ^2bju 2 - jh ,, (2.21) 

3=0 

where, naturally, u 2 -ih are the solutions to (|2.ip . with 2~^h in place of h, 

(&o,&i,...A) == (1,0,0,...,0)^- 1 (2.22) 
and V -1 is the inverse of the Vandermonde matrix with entries 
yij . = 2 -(i-l)(i-l) ) ij = l 5 ... )fc + L 
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The following result is a simple corollary of Theorem 12.31 

Theorem 2.5. In each situation when Theorem \2.3\ is applicable we have 
that the estimate 

\u h (t,x)-u^(t,x)\<N(\\f\\ m + \\g\\ m )h k+1 (2.23) 

holds for all (t,x) 6 Ht, h S (0, ho], where N is a constant depending only 
on d, m, |Ai|, Mq, . . . , M m , and T. 

Proof. By Theorem 12.31 

fe 



U 2 -ih = U 



(0) 



+ Ei tt(,)+f ^ j = o,i,...,fc, 



i!2* 

i=l 



with r 2 - Jft := 2 ^ k+1 ^r 2 -jf l , which gives 



«fc = E & ^ = £ & > (0) + EE *;^« w + E ^ 2 -^ fc+1 

3=0 j=0 3=0 i=l ' 3=0 

= - (0) + E £« w E^+E ^ = «w + E ft^^ 1 , 

i=l ' 3=0 3=0 3=0 

since 

fc fc 

^^ = 1, J^6j2~ y = 0, i = l,2, ...fc 

3=0 3=0 

by the definition of (bo, b^). Hence, 

k 

sup |u h , - «(°)| = sup | bjf 2 - jh \h k+1 < N(\\f\\ m + \\g\\ m )h k+1 , 
and the theorem is proved. □ 



Sometimes it suffices to combine fewer terms U2-jh to get accuracy of 
order fc + 1. To consider such a case for odd integers k > 1 define 

fe 

% = E > ( 2 - 24 ) 

3=0 

where 

(boM,...,^) := (1,0,0, ....OJV- 1 , fc = *fi, (2.25) 
and y _1 is the inverse of the Vandermonde matrix with entries 
yij :=4 -(i-i)(3-i) ) ij = 1,...^ + !. 

Theorem 2.6. Suppose that the assumptions of Theorem \2.3\ are satisfied 
and condition (Hi) is met. Then for Uh we have 

sup|«(°) -u h \ <N(\\f\\ m + \\g\\ m )h k+1 
for all h 6 (0, ho], where N depends only on d, m, M , . . . , M m , and T . 
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Proof. We obtain this result from Theorem 12.31 by a straightforward modi- 
fication of the proof of the previous result, taking into account that for odd 
j the terms with h? vanish in expansion (|2.18p when condition (iii) holds in 
Theorem 12.31 □ 

Example 2.1. Assume that in the situation of Theorem 12 .61 we have m = 8. 
Then 

Uh ■= |lifc/2 - \v>h 

satisfies 

sup|u (0) -u h \ < Nh A 

for all h e (0,ho\. 

The above results show that if the data in equation (|2.8|) are sufficiently 
smooth, then the order of accuracy in approximating the solution tS ' can 
be as high as we wish if we use suitable mixtures of finite difference approx- 
imations calculated along nested grids with different mesh-sizes. Assume 
now that we need to approximate not only but its derivative D a u^ for 
some multi-index a as well. What accuracy can we achieve? The answer is 
closely related to the question whether the expansion 

T j 

D a u h (t,x) = D a u {0) {t,x) + V — D a u^(t,x) + h k+1 D a r h (t,x) (2.26) 

i<?'<fc 

holds for all (t,x) € Ht and h £ (0, ho], such that 

\D a r h (t,x)\<N(\\f\\ m +\\g\\ m ) (2.27) 

for all (t,x) G H T , h £ (0, ho]. 

The result concerning this expansion and the following series of results 
appeared after the authors tried to extend the above theorems from the 
parabolic to the elliptic case. The main and rather hard obstacle is that the 
constants in our estimates depend on T and, actually, may grow exponen- 
tially in T. By the way, this obstacle is caused by possible degeneration of 
our equations and exists even if we consider equations in bounded smooth 
domain. 

To be able to give some conditions under which this does not happen, 
we introduce new notation and investigate smoothness properties of Uh with 
respect to x. As a simple byproduct of this investigation we also obtain 
smoothness of Uh with respect to h, which, by the way, cannot be derived 
from (EnSH - 

Take a function t\ defined on Ai taking values in [0, oo) and for A £ Ai 
introduce the operators 

T ht xv(x) = (p(x + /iA), 5 hj x = rxh' 1 ^^ - 1). 

Set 

l|Aill 2 =£MI 2 . 

AGAi 
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For uniformity of notation we also introduce A2 as the set of fixed distinct 
vectors i 1 , ...,£ d none of which is in Ai and define 

<W* = r oA, T hj£i = 1, A = Ai U A 2 , 

where To > is a fixed parameter. For A = (A 1 , A 2 ) € A 2 introduce the 
operators 

Th,\ = T h , x iT h \ 2 , 5 h: x = Sh,x^h,x 2 - 
For k = 1, 2, fj, € A k we set 

AgAi AgAi 

AgA AgAi 

Below B(R. d ) is the set of bounded Borel functions on M. d and .ft is the set 
of bounded operators K,\ t = K,h(t) mapping B(M. d ) into itself preserving the 
cone of nonnegative functions and satisfying /QT < 1. 

Finally, fix some constants 5 E (0, 1] and K £ [1, 00). 

Assumption 2.4. There exists a constant cq > such that c > c$. 

Remark 2.3. The above assumption is almost irrelevant if we only consider 
(|2.ip on a finite time interval. Indeed, if c is just bounded, say |c| < C = 
const, by introducing a new function v(t,x) = u(t,x)e~ 2Ct we will have an 
equation for v similar to (|2.ip with L9u — (c + 2C)v and fe~ 2Ct in place of 
L/jii and /, respectively. Now for the new c we have c + 2C > C. 

Assumption 2.5. We have m > 1 and for any /i € (0,/io], there exists an 
operator ICh = )Ch,m € .ft, such that 

< (1 - S) Qh($h,\<f) + KQ h (<p) + 2(1 - S)dC h ( l4,A^| 2 ) 
agA agA 

(2.28) 

on Ht for all smooth functions (p. 

Assumption 2.6. We have m > 2 and, for any /i € (0, ho] and re = 1, m, 
there exists an operator ICh = ICh,n € -ft) such that 

j^gA AgA 2 AgA 2 

+ J?i;^W) + 2(l-^(E |4,a^ 2 )+A-/C(^|4,a^ 2 ) (2.29) 

AGA AGA 2 AgA 

on fly for all smooth functions p. 

Obviously Assumptions 12.51 and 12.61 are satisfied if q\ and p\ are inde- 
pendent of x. In the general it is discussed in [5], the above as- 
sumptions impose not only analytical conditions, but they are related also 
to some structural conditions, which can somewhat easier be analized under 
the symmetry condition (S). 
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Assumption 2.7. For all t € [0,T] 

Xqx(t,x) is independent of x. (2.30) 

AeAi 

In the main case of applications we will require the last sum to be iden- 
tically zero as in Assumption [ 



Remark 2.4. Assumptions [231 and [2U1 are discussed at length and in many 
details in [5] and [6], and sufficient conditions, without involving test func- 
tions ip are given for these assumptions to be satisfied. In particular, it is 
shown in [6] that if condition (S) holds, m > 2, t\ = 1, Assumptions 12 . ll and 
12.21 are satisfied, and q\ > k for a constant k > 0, then both Assumptions 12 . 51 
and 12.61 are satisfied for any cq > and 5 € (0, 1), if ho is sufficiently small 
and To, K, and KL are chosen appropriately. Moreover, the condition k > 
can be dropped, provided, additionally, that cq is large enough (this time 
we need not assume that h is small). Remember, that by Remark 12.31 the 
condition that cq be large is, actually, harmless as long as we are concerned 
with equations on a finite time interval. Mixed situations, when c is large 
at those points where some of q\ can vanish are also considered in [6]. 

In [5] we have seen that Assumption 12.51 imposes certain nontrivial struc- 
tural conditions on q\ which cannot be guaranteed by the size of c$ if q\ is 
only once continuously differentiable. In contrast, even without condition 
(S), given that Assumptions 12. i l l2.5[l2~7H are satisfied and m > 2, as is shown 
in [6], Assumption 12.61 is also satisfied if cq is large enough. 

Theorem 2.7. Let Assumption \2. 1\ through \2. 6\ hold with m > 3. Let k > 
and I £ [0, m] be integers. Then for every multi-index a such that \a\ < I the 
function D a Uh is a continuous function on Ht and expansion (|2.26j) holds 
with D a rh satisfying (|2.27j) . provided one of the following conditions is met: 

(i) m>3k + 3 + l; 

(ii) m > 2k + 3 + I and condition (S) holds; 

(Hi) k is odd, m>2k + 2 + l, and conditions (S) and (I2.16P are satisfied. 
In each of the cases (i)-(iii) the constant N depends only on d, m, 5, K, tq, 
cq, |Aij, ||Ai||, Mo,... ,M m . In case (Hi) we have vk>> = for all odd j in 
the expansion. 

We prove this theorem in Section [71 Remember the definition of u~h and 
in (|2.2ip and (|2.24j) . The following is an obvious consequence of Theorem l2.71 



Corollary 2.8. Suppose that the assumptions of Theorem 2.1 are satisfied. 
Then 

sup \D a u h - W°>| < Nh k+ \\\f\\ m + \\g\\ m ), 

and if condition (Hi) is met then 

sup \D a u h - < Nh k+1 (\\f\\ m + \\g\\ m ), 

where N depends only on on d, m, 5, K, tq, cq, |Ai|, ||Ai||, Mq, . . . ,M m . 
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Remark 2.5. Observe that for k = Theorem 12.71 implies that 

sup \D a u h - D a u^ | < Nh (2.31) 

H T 

if m > 3 + \a\ and Assumption 12.11 through 12.61 hold. In addition one can 
replace D a Uh in (|2.3ip with 5%, where 

8? = S? 1 ■ ... • 5? d 

ft ft,ei h,£d 

and ej is the ith basis vector in M d . This follows easily from the mean 
value theorem and Theorem 12.91 b elow . The reader understands that similar 
assertion is true in case of Corollary 12.81 with the only difference that one 
needs larger m and better finite-difference approximations of D a . 

Next we investigate the smoothness of in x and h. Recall that for 
functions (p depending on h we use the notation D r h (p for the r-th derivative 
of ip in h. As usual, D®ip := ip. 

Remark 2.6. Suppose that Assumption 1 2 . 1 1 is satisfied. Take an hi € (0, ho), 
consider equation (12. ip as an equation about a function u ft (f, x) as function of 
(h,t,x) G [h\, ho] x Ht and look for solutions in the space 25 m (/ii) = 2$^ {hi) 
which is defined as the space of functions on [hi , ho] x Ht with finite norm 

sup \D a D r h u h (t,x)\. (2.32) 

\a\+3r<m l h l> h °]xH T 

It is obvious that the integrand in (|2.ip can be considered as the result 
of application of an operator, which is bounded in 23 m (/ii), to Uh(s,x). 
Therefore, a standard abstract theorem on solvability of ODEs in Banach 
spaces shows that there exists a solution of (|2,ip in Q3 m (/ii). Since just 
bounded solutions are uniquely defined by (|2.ip . we conclude that our Uh 
belongs to 25 m (/ii) for any hi € (0, ho)- Obviously, if the derivatives of the 
data are continuous in x, the same will hold for 

The above argument, actually, works if we replace \a\ + 3r < m with 
[a| + r < m in (I2.32p . We talk about ()2.32p in the above form because we 
will show that under our future assumptions the quantity (|2.32p is bounded 
independently of hi. 

Theorem 2.9. Let k > and m > 2 be integers and suppose that As- 
sumptions \2.1\ through \2.6\ are satisfied. Then, for each integer r > such 
that 

3k + r < m, 

the generalized derivatives D r D^Uh exist on (0, ho] x Ht, are bounded and 
we have 

\D r D k h u h \<N(\\f\\ m +\\g\\ m ), (2.33) 

where N is a constant depending only on m,6,co,To,K, Mo, M m , |Aj| ; 
and ||Ai||. In particular, G 53 m and 

ll^ft ||m m J ■ 
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We prove this theorem in Section El and in Section [6] we show that the 
following fact, used when we come to the elliptic case, is a simple corollary 
of it. 

Theorem 2.10. Suppose that Assumptions \2. I\ through \2. 61 hold with m > 
2. Then the constant N in (|2,10p depends only on m,S,CQ,To,K, Mq, 
M m , |Ai|, and ||Ai|| (thus, is independent ofT). The same is true for the 
constants N in Theorems W2 \ [Q j [Q j andVUk 

Additional information on the behavior of D r D^Uh for small h is provided 
by the following result which we prove in Section 

Theorem 2.11. Let k > 1 be an odd number and suppose that Assumptions 
\2.1\ through \2. 61 hold with m>2>k + l. Assume that the symmetry condition 
(S) and (|2.16p are satisfied. 

Then, for any integer r > such that 

3k + r < m — 1 

we have 

sup \D r D k h u h \ < N(\\f\\ m + \\g\\ m )h (2.34) 

for all h G (0,/to], where N depends only on m, 5, cq, tq, K, |Ai|, ||Ai||, 
M ,..., Af m . 

3. Main results for elliptic equations 

Here we assume that p\, q\, c, and / are independent of t and turn now 
our attention to the equations 

L h v h (x) + f{x) = xeR d , (3.1) 

+ f(x) = x € M d . (3.2) 
Naturally by a solution of (|3.2p we mean a function u on M rf such that 
it belongs to OS 2 and (|3.2p holds almost everywhere. Clearly, if a solution 
v belongs to 23 3 and q\, p\, c, and / are continuous functions on M. d , then 
(|3.2p holds everywhere. 

First we prove the existence and uniqueness of the solutions of equations 

and dS5D - 

Theorem 3.1. Suppose that Assumption \2. 1\ is satisfied with anm>0 and 
let Assumptions UTE and \2~^\ hold. Then equation (|3.ip /ias a unique bounded 
solution Vfi- Moreover, belongs to !B m . 

Proof. Observe that ()3.ip is equivalent to 

v h (x) = h 2 i(x)f{x) + £(x) ^2 X\Vh{x + Xh), 

AeAi 

where 

AeAi 
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It is seen that the existence and uniqueness of bounded solution of (|3.ip 
follows by contraction principle. Using smooth successive iterations yields 
that v h € <B m . □ 

Theorem 3.2. Let Assumptions \2.1\ through \2.bA hold with an m > 2, Then 
equation f)3.2f) has a unique solution v in the space OS 2 . Moreover, v £ QS m 
and i/iere is a constant N depending only on m, 5, cq, tq, K, Mq, M m , 
| Ai|, and ||Ai|| snc/i i/iaf 

IH|m<JV||/||m- (3-3) 

Proof. First we prove uniqueness. Let d € 23 2 satisfy (|3.2|) with / = 0. Take 
a constant ^ > 0, so small that c — f > co/2 and conditions (|2.28p and (|2,29p 
hold with c — v and <5/2 in place of c and 5, respectively. Then for each T > 
the function u(t, x) := e ut v(x)), (t, x) € Ht, is a solution of class 23^ of the 
equation 

d 

—u=(C + v)u on H T (3.4) 

with initial condition u(0,x) = v(x). Hence by virtue of Theorem 12.101 for 
every T > 

e uT \v{x)\ = \u(T,x)\ < N\\v\\ 2 , 

where N is independent of (T,x). Multiplying both sides of the above in- 
equality by e~ vT and letting T^oowe get v = 0, which proves uniqueness. 

To show the existence of a solution in <B m , let u be a function defined on 
Hoc such that for each T > its restriction onto Ht is the unique solution 
in 53™ of (|3.4p with initial condition u(0, x) = f(x) (see Theorem 12. ip . By 
Theorem [2710] 

sup^ \D r u\ <N\\f\\ m 

Ho ° r<m 

with a constant N depending only on m, 5, cq, to, K, Mq, M m , |Ai|, and 
||Ai||. Hence 



oo 



v(x) := / e u u(t,x)dt, xG 
J o 

is a well-defined function on v & 53 m , and 

/>oo 

Cv(x) = / e~ vt Lu(t,x)dt 



o 



1/t' a 



e vt (—u(t,x)-uu(t,x))dt = -f(x), 

where the last equality is obtained by integration by parts. Consequently, v 
is a solution of (13.41) and it satisfies estimate (13.31). □ 



Theorem 3.3. Let k > and suppose that Assumptions \2.1\ through \2.6\ 
are satisfied with an m > 3k. Then, for any h £ (0, ho] and for each integer 
r > 0, such that 

3k + r < m, 
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for the unique bounded solution Vh of (|3.1f) we have 

sup \D r D k v h \ < N\\f\\ m , (3.5) 

(0,h. ]xR d 

where N is a constant depending only on m,6,CQ,To,K, \Ai\, ||Ax||, Mq, 
M m . In particular, 

\\Vh\\m < N\\f\\ m . 



Proof. To prove (|3.5p , take a constant v > as in the proof of Theorem 13.21 
define u(t, x) := Vh(x)e vt , and observe that u is the unique bounded solution 
of 

d 

—u = L° h u - (c - v)u + e vt f, u(0, x) = v h (x). 

By Theorem 12.91 for any T > 

e vT \D r D k h v h {x)\ = \D r D k h u{T,x)\<Ne vT \\f\\ m + N\\v h \\ m , 

where N is a constant, depending only on m, <5, Co, To, K, |Ai|, ||Ai||, Mo, 
M m . By multiplying the extreme terms by e~ vT and letting T — > oo, we get 
the result. □ 

From estimate (|2.34p we obtain the corresponding estimate for the deriva- 
tives Of Vh- 

Theorem 3.4. Let the conditions of Theorem \2.11\ hold. Then for any 
integer r > such that 

3k + r < m — 1, 
for the solution Vh of ()3.ip we have 

swp\D r D k h v h \<N\\f\\ m h 

R d 

for all h £ (0,/io], where N depends only on m,5,co,TQ, K , |Ai|, ||Ai|| and 
Mo,..., M m . 



Proof. This theorem can be deduced from Theorem l2.11l in the same way as 
Theorem 13.31 is obtained from Theorem 12.91 □ 



Now we want to establish an expansion for v^, i.e., to show for an integer 
k > the existence of some functions v^°\...,v^ on and a function 
on R d for each h G (0, /i ] such that for all x £R d and /i G (0, ho] 

v h (x)=vW(x)+ Y ^ r v ( - j \x) + h k+1 R h (x), (3.6) 
i<j<fc J ' 

sup sup 1^1 < iV||/|| m (3.7) 
/iG(0,/i ] R d 

with a constant TV". 
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Theorem 3.5. Suppose that Assumptions \2. i\ through \2. 6\ are satisfied with 
an m > 3. Let k > be an integer. Then expansion (|3.6|) holds with 
being the unique Q3 m solution of (|3.2p and Rh satisfying (|3.7|) provided one 
of the following conditions is met: 

(i) m > 3k + 3; 

(ii) m > 2k + 3 and condition (S) holds; 

(Hi) k is odd, m > 2k + 2, and conditions (S) and (|2.16p are satisfied. 
In each of the cases (i)-(iii) the constant N in (|3.7p depends only on d, m, 
5,cq,tq, K , |Ai|, ||Ai||, Mo,...,M m . Moreover, when (Hi) holds we have 
yU) = o for all odd j. 

Proof. Take a small constant v > 0, as in the proof of Theorem 13.21 let u be 
a function defined on such that for each T > its restriction onto Ht 
is the unique solution in OS™ of 

d 

—u h = (L h + u)u h (t, x) 6 

u h (0,x)=f(x) x£R d , 
(see Remark \2.1\i . As in the proof of Theorem 13.21 we get that 

/■OO 

Vh(x) = / e~ ut u h (t,x)dt. 



By Theorem 12.31 in each of the cases (i)-(iii) we have 

u h (t,x)=u®(t,x)+ ^ j \t,x) + h k+1 r h (t,x), (3.8) 

for all (t, x) € H^, h € (0, ho], and by Theorem 12. 1UI we have 

k 

sup sup{K|+^|^')| + K|}<iV||/|U (3.9) 

he(0,ho] Hoc j_q 

with a constant N depending only on d, m, S,cq,tq,K, Mo,...,M m , \A±\ and 
||Ai||. Multiplying both sides of equation (|3.8p by e~ vt and then integrating 
them over [0, oo) with respect to dt, we get expansion (|3.6p with 

/•oo 

R h (x) := / e- ut r h (t,x)dt, 
Jo 

/•oo 

vW{x) := / e- ut u( j \t, x) dt, for j = 0, . . . , k. 
Jo 

Clearly, (|3.9p implies that (|3.7p holds with N depending only on d, m, 
5, cq, tq, K, Mo,...,M m , |Ai|, and ||Ai||. As we know the function in 
dUSD is the QS m solution of 

—u=(C + v)u (t,x)eH 00 , 

u(0,x)=f(x) x€R d , 
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which as we have seen in the proof of Theorem 13.21 guarantees that 

the unique 53 m solution of equation (|3.2p . □ 

Remark 3.1. We can show similarly that i = 1, ...,k, is the unique 
solution of the system 

3 

i=i 

in an appropriate class of functions (cf. Theorem 12.21) . 

The following result can be obtained easily from Theorem 12. 71 by inspect- 
ing the proof of the previous theorem. 

Theorem 3.6. Let p\, q\, c, and f satisfy the conditions of Theorem \3.5\ 
with m — l in place of m in each of the conditions (i)-(iii) for an integer I £ 
[0,m]. Then D a Vh is a bounded continuous function on M. d for every multi- 
index a, \a\ < I, and the expansion (|3.6j) is valid with D a Vh, {D a v^'}j =0 
andD a Rh in place ofvh, {v^}j=o and Rh, respectively. Furthermore, (|3.7p 
holds with D a Rh in place of R^ and a constant N depending only on d, m, 
5,cq,tq,K , |Ai|, ||Ai||, Mq, . . . , M m . In case (Hi) we have v^> = for all 
odd j in the expansion. 

Set 

k k 
j=0 j=0 

where (bo, b\, . . . , bk) and k, (bo, b\, . . . , b^) are defined in (|2.22p and in (I2.25|) . 
Then we have the following corollary. 

Corollary 3.7. Suppose that the assumptions of Theorem \3. 6\ are satisfied. 
Then for every multi-index a with \a\ < I, 

snp\D a v h -D a v^\<N\\f\\ m h k+1 , 

and if condition (Hi) is met then 

snp\D a v h -D a v^\<N\\f\\ m h k+1 , 

WL d 

where N depends only on on d, m, 5, K , tq, cq, |Ai|, ||Ai||, Mq, . . . , M m . 
4. Proof of uniqueness in Theorem 12. II and a stipulation 

We will see later that the proof of Theorem 12.31 only uses the existence 
of sufficiently smooth solutions of (|2.8p and (|2.13p . Therefore, if m > 3, 
uniqueness of u^ ' follows from expansion (|2.18p . If m = 2, one can use 
simple ideas based on integrating by parts. We briefly outline these ideas 
referring for details to [12], [13], |14] , 

First, one may assume that g = f = and let u^ ' be the corresponding 
solution. Then, by introducing a new function v = u(°)(cosh | | ) 1 one 
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reduces the issue to uniqueness of v, which satisfies an equation similar to 
(|2.5p with g = f = and different coefficients which we denote by q\, p\, 
and c = c, and, moreover, v,Dv,D 2 v € L2{Hj-)- After that one multiplies 
the equation for v by v and integrates over Ht- One uses integration by 
parts, and the fact that due to the assumption q\ > we have |-D<?a| 2 < 
4^asup \ D 2 q\\. One also uses Young's inequality implying that 

\v(d\q\)d x v\ < N\vql /2 d x v\ < q\{d\v) 2 + Nv 2 , 

and the fact that 2vp\d\v = p\d\(v) 2 . Then one quickly arrives at a relation 
like 

/ (N — c)\v\ 2 dxdt > / \v(T,x)\ 2 dx > 0, 
JH T J«. d 

where N is a constant independent of c. If c is large enough, the above 

inequality is only possible if v = 0, which proves uniqueness if c is large 

enough. In the general case it only remains to observe that the usual change 

of the unknown function taking v(t, x)e^ 1 in place of v for an appropriate A 

will lead to as large c as we like. 

Remark 4.1. Notice that apart from uniqueness in Theorems 12.11 and 12.21 
all our other assertions and assumptions are stable under applying molli- 
fications of the data with respect to x. For instance, take a nonnegative 
C € Co°(M a! ) with unit integral, for e > define C e (x) = s~ d ((x/s) and for 
locally summable tfj(x) use the notation 

^ =iP*( £ . 

Then q^f\p^\c^ , f^ e \ and will satisfy the same assumptions with the 
same constants as the original ones and will be infinitely differentiable in x. 

It is not hard to see that if our assertions are true for the mollified data, 
then they are also true for the original ones. For instance, let v e be the 
solution of (|2.5p with the new data. The uniform in e estimates of the 
derivatives in x and the equation itself, guaranteeing that the first derivatives 
in time are bounded, show that v e are uniformly continuous in [0, T] x {|x| < 
R} for any R. Then there is a sequence e n [ such that v Sn converges 
uniformly in [0, T] X {\x\ < R} for any R to a bounded continuous function v. 

This along with uniform boundedness of |Z) a u e |, |a| < m, lead to the fact 
that the generalized derivatives |-D a v|, |a| < m, are bounded and admit the 
same estimates as those of v e . Also since D a v £n — > D a v in the sense of 
distributions and all of them are uniformly bounded, we conclude that this 
convergence is true in the weak sense in any L2([0,T] x {|x| < R})- Now it 
is easy to pass to the limit in equation (|2.9p written for modified coefficients 
and v £ in place of u concluding that since the derivatives converge weakly 
and qy — > q\,..., — > / uniformly on Ht, v satisfies (12.9j) . 

Similar argument takes care of Theorem 12.21 (in which uniqueness will be 
derived from uniqueness in Theorem 12. ip . 

Our claim about stability of other results is almost obvious and 



ACCELERATED SCHEMES 19 

from this moment on we will assume that the data are as smooth in x as we 
like. 

5. Proof of Theorems 12.91 and 12.111 

In [5] (see there Theorems 2.3 and 2.1 and Corollary 3.2 if m = 0) and 
[6] we obtained the following result on the smoothness in x of the solution 
Uh to equation (|2.ip . 



Theorem 5.1. Suppose that Assumptions \2. II and \2.J\ are satisfied. Suppose 
that (i) if m = 1, then Assumptions UTM and \2.5\ are satisfied, and (ii) if 
m > 2, then Assumptions \2.5l \2.6l and 2.1 are satisfied. Then for 
h € (0, ho] we have that D k Uh, k = 0, ...,m, are continuous in x and 



SUpJ^ |D fc U h | < N i F m + G m ), (5.1) 

where 



Ht fc=o 



F n = J2 sn P\ Dk fh\, G n = ^ sup \D k g h \, 
k<n H ? k<n Rd 

and N depends only on To, m, 5, Co, K, |Ai| ; ||Ai||, Mo, M m (N depends 
on fewer parameters if m < 1). 

To proceed further we need a few formulas. 

Lemma 5.2. Let <p be a function on Ht and n > be an integer. 

(i) Assume that the derivatives of ip in x G M. d up to order n + 1 are 
continuous functions in x. Then for each h > 

Dh PA^AV = / 0nd \ +1 v(t, x + hOX) dO (5.2) 

AeAi AeAi Jo 

on Ht, where d\(p is introduced in (|2.12j) . 

(ii) Assume that the derivatives of ip in x up to order n + 2 are continuous 
functions in x, and that Assumption \2. 3\ holds. Then 

Dl h- l qx5h,\V = X) «A Al - 0)6 n d n x + \{t, x + MA) d8, (5.3) 
AeAi AeAi ^° 

on Ht. 

Proof. By Taylor's formula applied to ip(t,x + h8X) as a function of € [0, 1] 
^hx ( f(tj x )= / d\<p(t,x + h6\) d6 



Jo 

and 

-l 

52 



<5 M y(t, x) = d X ip(t, x) + h (i - e)a^(t, x + hex) do. 
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Multiplying the first equality by p\ and summing up in A over A\ we obtain 
(|5.2p for n = 0. Multiplying the second equality by q\, summing up in A 
over Ai we obtain (|5.3p for n = since 

q\d\v = o 

AGAi 

due to Assumption 12.31 

After that it only remains to differentiate n times in h both parts of the 
particular case of formulas (|5,2p and (|5.3p . The lemma is proved. □ 

Introduce 

= D{u h 

and observe that by Remark 12.61 under Assumption 12.11 the functions d™u J 
are well defined if n + j < m. By combining this with Lemma 15.21 and the 
Leibnitz formula we obtain the following. 

Corollary 5.3. Let Assumptions \2. il and \2.3\ be satisfied. Let k > 1 be an 
integer such that k + 2 < m. Then 

u[ k \t,x)= [ {L h u[ k \s,x) + R k h (s,x))ds (5.4) 
Jo 

on (0, ho] x Ht, where 

R k h (t,x) = Y,Cl Y, f O^px^x^ut^x + heX) 

i=l AgAi Jo 

+(i - 9) qx (t,x)(d x +2 u {k ~ i) )(t,x + hex)] dd. 

Now we are ready to prove Theorems 12.91 and 12.111 

Proof of Theorem \2.9l If m = 2 or k = 0, our assertion follow directly 
from Theorem 15. 11 Therefore, in the rest of the proof we assume that m > 3 
and k > 1. 

We will be using f)5.4f) . Observe that if 1 < i < /c, then 

(i + 2)+r+(fc-i)=& + 2 + r<3fc + r<m. 

Thus by Remark l2.6l we know that D %+2+r u^ ^ are bounded and continuous 
on ffy It follows that R k € 25 r . By Theorem 15.11 with r in place of m we 
obtain 

ife. := sup^ \Di U W\ < Nsu V J2\ DjR h\- 

Ht ■ ^ Ht ■ 

j<r 1 j<.r 

It is not hard to see that 

k i+2+j k 

H T j =1 n=1 i=1 
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Hence, 

k 

Ikr < Ny^^I k _ i: i +2+ r- 
i=l 

Here on the right the first index of I). r is reduced by at least 1 and the sum 
of indices increased by 2. Therefore, after k iterations we will come to the 
inequality 

Ikr < NIo,k+2k+r- 

It only remains to observe that Io,3fc+r < Io,m and the latter quantity is 
estimated in Theorem 15.11 The theorem is proved. 

Proof of Theorem \2.11\ First of all observe that the symmetry assumption 
and (|2,16p imply that for any smooth function <f(x), odd i > 0, and any 
multi-index a, such that \a\ < m, we have 

(D a Px)d\ +1 v = Yl (D<*q x )d{+ 2 V = 0. (5.5) 

AeAi AgAi 

If k = 1 and an integer n < r, then owing to (|5.5p 

\D n qx(t,x)(dlu h )(t,x + h8\)\ 
AeAi 

= \° n S ?A(t,i)[(^h)(M + WA)-^u h (t,i)]| 
AeAi 

< iV/tsup V < JV/i||u|| m < iV(||/|| m + \\g\\ m )h =: iVJ/i, 

where the last two estimates follow from the fact that r + 4 = r + 3/c + l <m 
and from Theorem 12.91 respectively. Similarly, 

\D n Px(t,x)(d 2 x u h )(t,x + h9X)\ 
AeAi 

= \D n Y,Px{t,x)[{dlu h ){t 1 x + heX)-dlu h {t 1 x)] \ <NJh. 
AeAi 

Hence, 

sup^|Z> n ^| < N(\\f\\ m + \\g\\ m )h < NJh 

^ T n<r 

and applying Theorem 15. II to (j5.4|) yields (|2.34[> . 

Now we proceed by induction on k. Assume that for an odd number 
j estimate (|2.34p holds whenever 3k + r < m — 1 and odd k < j. This 
hypothesis is justified by the above for j = 1 and to prove the theorem it 
suffices to show that the hypothesis also holds with j + 2 in place of j. Take 
an odd k and an integer r such that 

k <j + 2, 3k + r < m - 1 
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and again use (|5.4p . As above, to obtain (j2.34H it suffices to prove that 

supV \D n R k h \ < NJh. (5.6) 

1 n<r 

Take an integer n < r. Observe that if 1 < i < k and i is even, then k — i is 
odd and k — i<j + 2 — i<j and 

3{k-i)+i + 2 + n = 3k + n-2i + 2<m-l-2i + 2<m-l 

so that by the induction hypothesis 

sup ID™ V q x {t,x){d^ 2 u^~ i] ){t,x + h9\)\ < NJh. (5.7) 
H * AeAi 

If 1 < i < k and z is odd, then i + 2 is odd too and as in the beginning of 
the proof 

\D n qx(t,x)(di +2 ui k - l) )(t,x + h8\)\ 

AGAi 

= \D n Y qx{t,x)[{d^ 2 uf^t,x + he\)-df 2 u { t i \t,x)]\ 
AeAi 



<Nhsup Y \ D 



Ht i<k,Kr 



l+i+3,,{k-i)\ 



where the last sup is majorated by NJ owing to Theorem 12.91 since 

3(k — i)+r + i + 3<m — 1 — 2i + 3 < m. 

In both situations we have (|5.7p . Similarly, if 1 < i < k and i is odd, then 
i + 1 is even and 

\D n Y PA(t,x)(^ +1 4 fc - i) )(t,x + ^A)| 
agAj 

= \D n Y Px(t,x)[(d^ 1 u ( £- i) )(t,x + hOX) - dl +1 u ( *- l \t,x)] | 

AGAi 



<Nhsup Y \ D 



Ht i<k,Kr 



l+i+2,Xk-i)t 



where the last sup is majorated by NJ again owing to Theorem 12.91 since 

3{k-i)+r + i + 2<m-l-2i + 2<m. 

Finally, if 1 < i < k and i is even, then k — i is odd, k — i < j + 2 — i < j, 
and 

3(fc-i)+r + i + l<m-l-2i + l<m-l, 
so that by the induction hypothesis 

\D n Y Px(t,x)(d i x +1 u (k ~ i) )(t,x + h9X)\ < NJh, 

AGAi 

which is now shown to hold in both subcases. By combining this with (|5.7p 
we come to (15.6P and the theorem is proved. 
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6. Proof of Theorems EH E21 and LTTOl 

Proof of Theorem \2.1\ First we replace q\ with symmetric ones using the 
fact that the symmetrization does not affect formula (12. 7ft . To this end 
introduce 

A? = Ain(-Ai), Ai = AiU(-Ai). 

On Af we set q± x = (1/2)( 9a + <Z-a)- If A G ±(Ai \Af ) we set q x = (l/2)g ±A . 
Then Ai and q\ satisfy the symmetry condition (S) and can be used to 
represent the first term on the right in (|2.7p in place of the original ones. 
Next, we redefine and extend p\ introducing p\ on Ai, so that p\ = Mq+p\ 
on Af , for A € Ai \ Af we set p±\ = Mq ± {l/2)p\, and for —A G Ai \ Af 
we set p±\ = Mq =p (l/2)p_ A . (Remember that for the constant Mq from 
Assumption 12.11 we have \p\\ < Mq.) Then Ai and p\ can be used to 
represent the second term on the right in (|2.7p in place of the original ones. 
One of the advantages of the new p\ is that p\ > 0, which implies that the 
new x\ satisfies Assumption 12.21 

Define t\ > arbitrarily. As in Remark 6.4 of [5] and Remark 4.3 of [6] 
one shows that Assumptions 12.51 and 12.61 are also satisfied for any 5 G (0, 1), 
say 5 = 1/2, if c is sufficiently large (independently of h) and To > 0, K, and 
fC are chosen appropriately and depending only on d, | Ax [, ||Ai ||, Mq, Mi, M2. 
We first concentrate on the case that c is indeed sufficiently large. In that 
case by Theorem 12 .9\ for h G (0, ho], there exists a unique solution Uh(t,x) 
of class 23^ satisfying equation (|2.1h with Lh in place of L^, where Lh is 
constructed from Ai, q\, and p\. Furthermore, 

IMU<iV(||/|| m + ||2|| m ), (6.1) 

where N is a constant depending only on m,infc, |Ai|, Mq,..., M m , and 
||Ai||. Upon observing that owing to Remark 12.21 

\L h u h \ < N(sup \D 2 u h \ + sup \Du h \ + sup 1^1) 

H'j' Hj' Hj' 

with independent of h, we conclude from the equation for that their 
first derivatives in t are bounded uniformly in h. Therefore, there exists a 
sequence h(n) j such that u^ua converges uniformly on [0, T] x {x : \x\ < 
R} for any R to a continuous function v. Then (16. If) implies that v G 53 m 
and 

IM| m <iV(||/|| m + ||<7|| m ) (6.2) 

with the same N as in (16. ip . If we take t a = 1, then Remark 6.4 of [5] and 
Remark 4.3 of [6] imply that both N's can be chosen to depend only on d, 
m, inf c, |Ai|, and Mo, M m . 

Next, the modified equation (|2.ip yields that for any <p G Co°(R d ) and 
t G [0,21 

/ Uh(t,x)(j)(x) dx = / g(x)(j)(x)dx 

jR d jR d 
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+ / / u h(s,x)[(l/2)A hX (q\(l)) +8 h x{Px4>)\{s,x)dxds 

Jo J» d ri 

+ / / (—cuh. + f)4>(s,x)dxds. 

JO JR d 

We pass to the limit in this equation and find that v satisfies an integral 
equation, integrating by parts in which proves that v is a solution of ()2.8|) . 

Finally, we notice that the case that c is not large is reduced to the above 
one by usual change of the unknown function taking v(t, x)e xt in place of v 
for an appropriate A, which leads to subtracting Xv from the right-hand side 
of (|2.5p . For the new equation we then find a solution admitting estimate 
(|6.2p with N independent of T but coming back to the solution of the original 
equation will bring an exponential factor depending on T. 

This and uniqueness proved in Section [4] finish proving the theorem. □ 

Remark 6.1. In the above proof we considered arbitrary t\ > for the 
following reason. If Assumptions 12.11 through 12.61 hold with m > 2, then 
by Theorem 12.91 estimate (|6,ip and hence (I6.2p hold with iV depending only 
on m, 6, cq, To, K, Mo, M m , |Ai|, and ||Ai||. This proves the assertion of 
Theorem 12.101 regarding the constant N in Theorem 12.11 

Proof of Theorem \2.2l Notice that for each j = 1, . . . , k equation (|2.13p does 
not involve the unknown functions ttW with indices I > j. Therefore we can 
solve (|2.13p and prove the statements (i) and (ii) recursively on j. 

First we prove that there is at most one solution (w 1 ', . . . ,u^) in the 
space OS 2 x • • • x *B 2 . Denote 

3 

8=1 

We may assume that = 0. Then clearly S\ = and by Theorem 12.11 
we have = 0. If for a j € {2, . . . k} we have = = ■ ■ ■ = 
u^~^' = 0, then clearly Sj = which by Theorem 12.11 yields v&> = 0. Hence 
the statements on uniqueness follow because for every j = 1,2, ...,k we 
obviously have <B m_3j ' C *B 2 when m > 3k + 2 and W n ~ 2j C *B 2 when 
m > 2k + 2. 

While dealing with the existence of a solution first take j = 1. Observe 
that by Theorem 12.11 we have vf® € 25 m with m > 5 in case (i) and with 
m > 4 in case (ii). Thus in case (i) we have S± £ <B m ~ 3 C 25 2 and by 
Theorem 12.11 it follows that there exists € 5S m_3 satisfying (|2. 13[) and 
admitting the estimate 

||« (1) ||m-3<JV||u (0) || ro . 

Taking the estimate of the last term again from Theorem 12.11 we obtain 
(I2.14p for j = 1. In case (ii) we have actually better smoothness of S\, 
because the first sum in (|2.1ip is zero for % = 1 and, for that matter, for 
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all odd i. It follows that Si G <B m ~ 2 and this leads to (12TT5T) for j = 1 as 
above. By adding that under the conditions (S) and (12.16P we have = 0, 
Si = 0, and = 0, we obtain pTTl) for j = 1. 

Passing to higher j we assume that k > 2. Suppose that, for a j 6 
{2, we have found i/ 1 ),...,^ -1 -* with the asserted properties. Then in 

the case (i) we have 

£® u U-i) G <B m " 3j C <B 2 

for i = 1, . . . , j, since 

m - 3(j - i) - (i + 2) = m - 3j + 2i - 2 > m - 3j > 2. 

Hence 5j G «B m ~ 3 -? and therefore by Theorem 12.11 there exists «w) g <8 m -3j 
satisfying (|2.13p and admitting the estimate 

lk (i) IU- 3j <JvJ^||ttO--*)|| ro _ 3i+a 

i=l 

j 

< llu^lU-g^-i) < N(\\f\\ m + \\g\\ m ), 

i=l 

where the last inequality follows by the induction hypothesis. 

In case (ii) we take into account that due to condition (S) we have 

]T qx d\ +2 ip = 0, (6.3) 
AeAi 

and due to condition (|2.16p we have 

£ paS^V = (6.4) 
AGAi 

for odd numbers i and sufficiently smooth functions (p. It follows that in 
case (ii) for i = 1, j we have 

since ^WuCj-i) e $sm-2(j-l)-2 and for ^ > 2 

m - 2(j - i) - (i + 2) = m - 2j + % - 2 > m - 2j > 2. 

Hence Sj G <B m_2j and therefore by Theorem [2TT1 there exists G <B' m - 2 -? 
satisfying (|2.13p and admitting the estimate 

3 

\m-2j+3 

i=2 

v 11. Ci-ilii 

Im— 2(7— t)i 

i=l 

and by using the induction hypothesis we come to (|2.15p . 

Furthermore, in case (ii) if (|2.16p is satisfied, our induction hypothesis 
says that «W = for all odd I < j - 1. If j is even, then, obviously, = 



\n^\\ m -2j < N\\u^\\ m - 2j+2 + 



<N^2\\u^~ 
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for all odd I < j as well. If j is odd then to carry the induction forward it 
only remains to prove that u^) = 0. However, for odd i we have 

£(o u (j-o = o 

due to (|6.3p - (|6.4p . This equality also holds if i > 2 and i is even, since then 
j — i is odd and u^~^ = by assumption. Thus, Sj = and = 0. □ 

Remark 6.2. The above proof is based on Theorem 12.11 and leads to esti- 
mates (|2.14p and (|2. 15|) with N depending only on the same parameters 
as in Theorem 12.11 Therefore, according to Remark 16.11 if Assumptions 12.11 
through 12.61 are satisfied and the restrictions on m and k from Theorem 12.21 
are met, then the constants iV in estimates (j2.14h and (|2.15D depend only 
on m, 6, co, to, K, Mo, M m , |Ai|, and ||Ai||. This proves the part of asser- 
tions of Theorem 12.101 concerning Theorem 12.21 The proof of its remaining 
assertions can be obtained in the same way and is left to the reader. 

7. Proof of Theorem [273] and [2~71 

We need some lemmas. The first one is a simple lemma from undergrad- 
uate calculus on Taylor's expansion. 

Lemma 7.1. Let F be a real-valued function on (0, 1] such that for an 
integer m > the derivative F^ m+1 \h) of order m+1 exists for all h € (0, 1], 
and F( m+1 ) is a bounded function on (0, 1]. Then 

FW(0) := limits) 

exist for < k < m, and 
holds for h € [0, 1] with 

Rm(h)= [ h ^-^F( m + 1 Hs)ds, 

Jo rn\ 

so that 

i. m+l 

\Rm{h)\ < sup \F^ m+l \s)\- ( — - for all h € [0, 1]. 

ae(o,i] (m + 1)! 

To formulate our next lemma we recall the operators Lh, C and 
defined in (|2.2p . (12. 7{\ . and (12.111) . respectively, and for each h S (0, ho] and 
integer j > introduce the operator 

0f=L h -C- 

l<i<j 
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Lemma 7.2. Let Assumption \2.3\ hold. Assume that for some integer I > 
the functions p\, q\ belong to for all X E A\. Then for any integer j > 

< iV|Mlz +J+3 ^ +1 (7-1) 

for all h € (0, /i ] and if G <B' +J+3 , w/iere N is a constant depending only 
on |Ai|,Mo,...,M,. 

Proof. We may assume that the derivatives in x of up to order / + j + 3 
are bounded continuous functions on Ht- By Lemma 15, 21 the derivatives of 
the function L^ip in h up to the (I + j + l)st order are bounded functions 
on (0, ho] x Ht and 

(£(f>)(t,x) = lim(L h <p)(t,x), 

(£®<l>)(t,x) = lim(DlL h 0)(t,x). 

Thus applying Lemma 17.11 to F{h) := Lhtp(t, x) for fixed (t, x) and using 
Lemma 15.21 we have 

= V 9a f ^—11. [\i-e)6i +1 &> x +3 p(t,x + tf6\)d6dtf 

AeAi ^ J ' ! ^ 
Now estimate (IT. If) follows easily. □ 

The next lemma is a version of the maximum principle for d/dt — Lh- It 
is a special case of Corollary 3.2 in [5]. 

Lemma 7.3. Let Assumption \2.1\ with m = be satisfied and let Xh,x > 
for all A £ Aj. Let v be a bounded function on Ht, such that the partial de- 
rivative dv(t, x)/dt exists in Ht- Let F be a nonnegative integrable function 
on [0, T], and let C be a nonnegative bounded function on Ht such that 

v := sup(C - c) < 0. 

Assume that for all (t, x) 6 Ht we have 

< L h v + Cv+ + F, (7.2) 

where v{t) = sup{t>(i, x) : x £ R d }. Then in [0, T] we have 

v(t) < v + (0) + \v[ 

where a + : = (ja| + a)/2 for real numbers a 



v(t) < «+(0) + |i/| _1 supF, (7.3) 

[o,t] 
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Proof of Theorem \2.3l By taking Uhe~^ M ° +1 " )t in place of Uh, we may assume 
that c > 1. Consider first the case k = 0. Since m > 3, by Theorem 12.11 
equation f|2.Tj) has a solution J ', which belongs to <B m and estimate (|2.10p 
holds. Clearly, w := u^ — vf^ is the unique bounded solution of the equation 
ft 

w(t, x) 



[ (L h w(s,x)+F{s,x))ds, (t,x)EH T , (7.4) 
Jo 

where F := of ] u^ = L h u^ - CvW . By Lemma E2 and estimate fl2~T0l) 
|rfV°)|| < TV ]T (IIpaIIo + hxh)\\u^hh < N(\\f\\ 3 + \\g\\ 3 )h 



'h 
,(0)| |(j 

AGAi 

with constants N depending only on d, \A±\ Mq, Mi, M3, and T. After 
that an application of Lemma 17.31 to equation (|7.4p proves the statement of 
Theorem 12.31 for k = 0. 

Let A; > 1 . Then by Theorem 12.21 the system of equations (|2.13p has a 
bounded solution {u^}^ =l . Observe that for 

k 1 j 

u; := it/,, — u^— (7-5) 
j=o J - 



we have equation (|7.4p with 



F := L^W - £«C°> + ^L h n(^ - - G, 



7 7 
i=i J i=i J 



and 

G 



k j -. k k 



E E MT^T^ti^V = E E TTTT- — —C^u^^^h? 
^ ^ - 1)\ ^ ^ - 1]\ 

j=l 1=1 yJ > i=l j=i yJ ' 



i=l i=0 ' i=0 t=l 



(0 U (0 



= Ett E ^ W « W - 

Hence by simple arithmetics 



7! ^— ' i\ 

j=0 J \<i<k-j 



p =ET°k N)tt0) - (7 - 6) 

Notice that 

A; — j + 3 < m — 3 j for j = 0, 1, . . . , k in case (i) , 
k — j + 3 < m — 2j for j = 0, 1, . . . , k in case (ii) , 
k — j + 3 < m — 2 j for j = 0, 1, . . . , k — 1 in case (hi). 
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Therefore by Theorem 12.21 under each of (i) , (ii) , and (hi) 

\\u^\\ k _ J+3 <N(\\f\\ m + \\g\\ m ) 
for j = 0, 1 . . . , k (u^ = in the case (hi)). Thus by Lemma 17.21 

\\O { h h - j) u^\\ < Nh k -J +1 \\u®\\ k - j+3 < Nh k+1 -i(\\f\\ m + \\g\\ m ). 
Consequently, 

ll^llo < N(\\f\\ m + ||9l|m)^ +1 for h G (0, ho], 

where A?" depends only on d, m, |Ai|, Mq, . . . ,M m , and T. Hence we get 
(|2.18p by Lemma [7T3l and the proof is complete. □ 

Proof of Theorem 2.7 Coming back to the above proof of Theorem 12.31 we 
see that function (|7.5p satisfies (|7.4p with F given by (|7.6p . We notice that 

k — j + 3 + I < m — 3j for j = 0, 1, . . . , k in case (i) , 
k — j + 3 + I < m — 2j for j = 0, 1, . . . , k in case (ii) , 

k — j + 3 + I < m — 2j for j = 0, 1, . . . , k — 1 in case (hi). 

Therefore by Theorem 12. 1\ when k = 0, and by Theorem 12 .2^ when k > 1, 
under each of (i), (ii), and (hi) 

||u^|| fc -i43+J <N(\\f\\ m + \\g\\ m ) 

for j = 0, 1 . . . , k (u (fc) = in case (hi)). By Theorem 12.101 the constant N 
depends only on m, 5, cq, tq, K, Mq, M m , |Ai|, and ||Ai||. By Lemma 

m 

where N is a constant depending only on |Ai|, Mq,. . . M\. Hence 

< JV(||/||m + \\9\\m)h k+1 for h G (0, h }. 

Consequently, applying Theorem 12.91 to equation (|7.4p , for any multi- index 
a, \a\ < I, for 

k 

we have 



7 

3=0 J 



r (»)|| = / l -(^l)|| J D- U ;|| <iV(||/||, 



\9\ 



with a constant N depending only on m, d, 5, Co, To, K, Mq, M m , |Ai| 
and ||Ai||, which proves the theorem. □ 
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